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We investigate the problem of propagation of three-component resonant hght pulses with adiabati- 
cally varying amplitudes through a medium consisting of atoms with the tripod level configuration. 
By means of both analytic and numerical methods we found the two modes of shape-preserving 
pulse propagation. The pulse propagation velocity is found to be either equal to the speed of light 
or significantly slowed down, depending on a particular propagation mode. 



PACS number: 42.50.Gy 

I. INTRODUCTION 

Coherent population trapping (CPT) is a well-known 
phenomenon of preparation of atoms in a coherent su- 
perposition of ground or metastable state sublevels (so- 
called dark state), which is immune to excitation by a 
two-component laser radiation under the two-photon res- 
onance condition [1]. Since the laser radiation is not 
scattered by atoms in the dark state, the radiation ab- 
sorption is dramatically reduced. This effect is called 
electromagnetically induced transparency (EIT) and is 
actively studied since early 90's [2]. One of the most 
striking features of EIT is possibility of shape-preserving 
propagation of light pulses with slowly (adiabatically) 
varying amplitudes at the group velocity significantly re- 
duced with respect to the speed of light in vacuum, c 
[3]. Shape-preserving electromagnetic pulses propagat- 
ing in a coherent atomic medium at the reduced group 
velocity were called in Ref. [3] "adiabatons" . Experi- 
mental observations of light a pulse group velocity less 
by many orders of magnitude than c has been repeat- 
edly reported [4]. Slowing down the laser light followed 
by spatial compression of the pulses provides a unique 
possibility for design of nonlinear-optical devices operat- 
ing on a few-photon level [5] . Extreme sensitivity of CPT 
and EIT to deviations from the two-photon resonance al- 
lowed to observe experimentally large Kerr nonlinearity 
[6] and absorptive optical switching [7] in cold rubidium 
vapor. Such nonlinear optical phenomena, along with 
the possibility of reversible conversion of a photonic ex- 
citation to a collective spin excitation [8] and trapping 
light in a medium with the photonic band gap induced 
by a periodic modulation of the EIT resonance [9] , are of 
great importance for quantum information storage and 
processing. 

A novel direction in CPT and EIT studies is related to 
the systems admitting more than one dark state for the 
given real Rabi frequencies, and phases, Xj, associ- 
ated with the resonantly driven transitions. The simplest 
scheme of such a kind is the tripod scheme displayed in 
Fig. 1. Stimulated Raman adiabatic passage (STIRAP) 



in such an optically thin atomic medium with the tri- 
pod level scheme was investigated theoretically [10] and 
demonstrated experimentally [11] by Bergmann and co- 
workers. A proposal to use the tripod scheme as a phys- 
ical implementation of a qubit has been made recently 
[12]. A 5-level scheme being the extension of the tripod 
scheme was considered in Ref. [13]. 
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FIG. 1. Tripod scheme of levels driven by resonant electro- 
magnetic fields. 

The very specifics of the tripod scheme is that during 
adiabatically slow change of the external field parameters 
transitions between the two dark state occurs. These 
transitions are described by a non-Abelian phase ma- 
trix [14], which is a generalization of a geometric (Berry) 
phase [15] to the case of degenerate eigenstates of an 
adiabatic Hamiltonian. One may expect that these tran- 
sitions give rise to rich and complicated dynamics of laser 
pulse propagation in an optically dense medium with the 
tripod level configuration. However, only few theoretical 
works on EIT in such media are available. Paspalakis and 
Knight [16] considered parametric frequency generation 
for the case of time-independent fields at the medium en- 
trance and calculated the group velocity of a weak probe 
field. Petrosyan and Malakyan [17] investigated theo- 
retically EIT in a tripod medium as a tool for optical 
cross-phase modulation and high-precision magnetome- 
try in the weak probe field limit. The value of the group 
velocity obtained in [16,17] is strongly reduced with re- 
spect to c in the same way as in the standard case of 
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EIT in a three- level medium [3]. In the theoretical in- 
terpretation of the experiment on four-wave mixing in a 
solid-state system with the tripod level configuration [18] 
and other numerical calculations by Ham [19] related to 
that system, small optical density of the medium was 
assumed. 

The aim of the present paper is to study pulse propa- 
gation in a medium with the tripod level scheme (here- 
after briefly called "tripod medium") for a general case, 
in which none of the three resonant electromagnetic fields 
is assumed to be weak compared to others. The paper is 
organized as follows. In Sec. II we present the set of equa- 
tions treating pulse propagation in the adiabatic regime 
in a tripod medium beyond the weak probe approxima- 
tion . In Sec. Ill the analytic solutions describing slow 
and fast pulse propagation are obtained. Sec. IV contains 
the rcsiilts of our numerical calculations and their inter- 
pretation. Sec. V deals with some particular regimes of 
propagations. Sec. VI is devoted to conclusive remarks. 



II. BASIC EQUATIONS 

If the three electromagnetic fields are tuned exactly in 
resonance with the corresponding transitions |j) <-> |0), 
j = 1,2,3, the Hamiltonian in the interaction represen- 
tation reads as 



^ = -n^fi,|0>(i|+H.c., 



(1) 



Oj = jfijje*'^' = dojEj/h, where doj is the dipole moment 
matrix element of the given transition. The electric field 
in the jth laser wave is Ej ex.p[ikj{z — t/cj] + c.c, kj be- 
ing the radiation wave number. The complex amplitude 
Ej is a slowly varying function of z and t. Expanding 

the atomic wave function as = ao |0) -|- Xlj=i % b)> 
we obtain the Schrodinger equation for the probability 
amplitudes: 



icij = - fl*ao, j = 1,2,3, 



iao = — ^jdj- 



(2) 



The set of shortened Maxwell equations for slowly vary- 
ing field amplitudes can be written as 



d_ ld_ 

dz c dt 



= iGjUoa*, j = 1, 2, 3, (3) 



where Gj = 27TkjndQj/h and n is the atomic number 
density. Taking into consideration propagation effects 
described by Eq. (3) is the essence of the theory devel- 
oped in the present Section, in contrast to the theory of 
Ref. [10] , which applies to the case of a refractively thin 
medium. 



Hereafter we assume that all the matter-field coupling 

constants are equal: 



Gi = G2 = Gs = G. 



(4) 



Violation of this assumption leads to adiabaticity break- 
down during the pulse propagation and subsequent pulse 
front steepening [20] . Thermal motion of atoms leads, be- 
sides reduction of the effective number density of atoms 
in resonance with the laser radiation, to a similar effect 
of pulse front steepening [21]. However, the pulse shape 
distortion effects manifest themselves at propagation dis- 
tances much larger than the typical propagation distance 
associated with an adiabaton-likc pulse formation [20,21]. 
Therefore we can neglect both the differences of oscilla- 
tor strengths associated with the three laser-driven tran- 
sitions and thermal motion of atoms. We also neglect the 
radiative decay of the excited state |0) since it plays no 
role in the adiabatic regime, because of negligible popu- 
lation of the excited state [3] . 

We parameterize the Rabi frequencies by introducing 

/ 3 \l/2 

the generalized Rabi frequency 57=1 J2j=i 1% 1^ ) 
two angular variables ip and 1?: 



01 = sini?cos(/5e'^^n, 

02 = cosiJcosy e'-^^fi, 
^3 = smcpe'^^n. 



(5) 



There are two mutually orthogonal non-absorbing (dark) 
states associated with the Hamiltonian of Eq. (1): 



$(1)^ = cos'&e-'^' |1) - siwde-'^^ |2) , 

$(2) ^ = sin ^smipe-^^^\l)+ (6) 
cosiJsini/je"'-^^ |2) — cos</?e~'-^^ |3) . 



An atom initially prepared in either of these two states 
remains unexcited since 



H = 



1,2. 



(7) 



An atom also remains unexcited if the parameters of the 
laser radiation vary in time slowly enough to satisfy the 
adiabaticity conditions 



and 



Xj^fl, i = 1,2,3. 



(8) 



(9) 



However, in the latter case there are transitions [14] be- 
tween the dark states defined by Eqs. (6), where the in- 
stantaneous values of the varying angles (f and 1? and 
phases Xj enter. If at t ^ —00 an atom was in the sth 
dark state, its wave function at subsequent instants of 
time is 
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s' = l 

where the matrix B obeys the equation 

B,,, (t) + ^ Ess" {t)As',„ (t) = 

s" 

with the initial condition 

Bss'( — OO) = ^ss', 

and 



_9 



(10) 

(11) 

(12) 
(13) 



ExpHcitly, 

^11 = -i{Xi cos^ '& + X2 sin^ i?), 

Ai2 = 1? sin — i{xi — X2) sin ■!? cos "& sin y>, 

^21 = 

A22 = -«[(xi sin^ ^? + X2 cos^ t?) sin^ + 



X3COS If]. 



(14) 



It is easy to show that if the phases of the laser fields 
are kept constant at the medium entrance, then 



Xi=0, j = l,2,3, 



(15) 



in the whole tripod medium. The opposite is not true. 
If the absolute values of the field amplitudes are con- 
stant at the medium entrance, but the phases are mod- 
ulated, then the absolute values of the fields amplitudes 
and, hence, 1? and become time-dependent inside the 
medium. In the present paper we consider only the case 
when Eq. (15) holds. In this case Eqs. (14) are reduced 
to 



An =^22=0 



Aio = -Aoi = V 



where 



V = if sm (/J 



(16) 



(17) 



and v{z, —00) = 0. Then Eqs. (11, 12) yield the following 
result [10]: 



-Bii = B22 = cosv, 



B12 = -B21 = sin u. 



(18) 



We assume that the tripod medium occupies the half- 
space z > 0. Initially, at t ^ —00, all the atoms in the 
medium are in the coherent superposition of the dark 
states 



|^(— 00)) = cos/3 



) -hsin/3 



$(2) 



(19) 



The boundary conditions for the fields at the medium 
entrance 0(0, f) = ^oit), 'd{0,t) = 'i)o{t), and f{0,t) = 
ipo{t) are consistent with Eq. (8). Thus the adiabatic 



regime of the laser radiation propagation inside the 
medium is ensured. It is convenient to introduce new 
variables ( = z and r = t — z/c, as in Ref. [3]. Re- 
spectively, the derivatives over the new variables are 
d/{dT) = d/{dt) and d/{dC) = d/{dz) + c-^d/{dt). 

Now we can solve self-consistcntly the set of 
Schrodinger — Maxwell equations (2, 3). First of all, 
we note that in the adiabatic regime ao is very small, 
and the probability amplitudes of the low-energy states 
{j = 1,2,3) are, according to Eqs. (10, 18, 19), 

aj = cos(!/ + p) {j\ $(1) ) + sin(!/ + (3) {j\ $(2) ) . (20) 

Then we find c;asily, that, similarly to the case of adia- 
batic pulse propagation in a A- medium [3] , 



(21) 



i.e., = ^o{t). Then we use the trick first applied in 
Ref. [3]: We express the small probability amplitude of 

the excited state as ao = —{i/Vl*)daj/{dT) and substi- 
tute this expression into the shortened Maxwell equa- 
tions (3). We get (9|f}j/0|V(9C) = {G/n'^)d\aj\^ /{Or), 
j = 1, 2, 3, or, explicitly, 

d d 

— (sin 1? cos ^pf = TT- fcos(i' -|- /?) cos 1? -|- 

oQ ow 

sin(z/ + (3) sin i9 sin ip\^ , 
d d 

— (cos ^ cos <p)^ = TT- \^ cos(i/ + 3) sin 1} + 

ow 

sin(:/ + 0) cos 1? sin tp]"^, 
^ sin^ ^ = ^[sin(zv + /?) cos^]'. (22) 

Here we introduced, instead of r, a new variable (nonlin- 
ear time) 



1 f^nl{r')dT', 



(23) 



which has the dimension of length. Then Eq. (17) takes 
the form 



dv d-d 
dw dw 



sm (f. 



(24) 



All the initial conditions set at t ^ —00 apply now to 
w; = 0. 

Only two of Eqs. (22) are independent. After some 
tedious calculations they are reduced to 

sin(t/ + P) /d_ _d_\ 
cos ip \d(^ dw J ^ 

cos(. + /?)(! + A). = 0, 

£^£(!i±^^+sin(. + /3)f =0. (25) 
cos ip oQ oQ 
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It is convenient now to change the variables to ui = 
( — w and U2 = w. The set of Eqs. (24, 25) takes the 
form 

cos(z^ + /3) dip 
cos (fi dui 
sin(i/ + (3) dip 



COS(p 

dv 



dU2 

dv 

dU2 



siii{iy + P)- — 

oui 


= 0, 


(26) 


^^ 9^ 
COS(!/ + /3)-— 

dU2 


= 0, 


(27) 


/d-d d'&\ 
\ dui du2 J 


shup. 


(28) 



III. 



SLOW AND FAST PULSES: 
SOLUTION 



THE ANALYTIC 



The set of Eqs. (26 - 28) is especially convenient for 

searching analytic solutions in a case when the unknown 
functions ip and i? depend on only one of the variables 
Wi, U2- We find two classes of solutions. 

The first one is the class of slow pulses. In this case 
the unknown functions depend only on ui = — w. 
The group velocity Vg of pulses of such has the same 
form as that of adiabatons in a A-medium [3]: Vg = 

{c-^+G/Q?y and can be much less than c. All the 
derivatives over U2 vanish, thus making Eq. (27) an iden- 
tity. The two remaining equations (26) and (28) become 
ordinary differential equations, yielding the general solu- 
tion in the parametric form: 



I cos 1^1 



Ci 



|cos(!/ + /3)p 



\sm{d-C2)\ 



\s\n{v + l3)\ 



(29) 



Here Ci, C2 arc arbitrary constants, and p is any func- 
tion of ^ — w compatible with the adiabaticity conditions 
(8). 

Similarly, wc find a general solution for the class of fast 
pulses, propagating at the speed of light: 



I cosi^l = 



C3 



sin(z/ + /3)|- 



|sin(i?-C4)| 



|cos(i^ + /3)| 



(30) 



Here C3, C4 are arbitrary constants, and p = piw) must 
be compatible with Eq. (8). 

Although the set of Eqs. (26 - 28) looks rather simple 
and symmetric, our attempts to find its general solution 
in the case of dependence of and •& on both u\ and U2 
have been unsuccessful. However, we can prove that a 
time-dependent solution in the parametric form 



ip = (p{p), 
does not exist if 

OUl 



'd = -dip) 



0U2 



(31) 



(32) 



Indeed, Eqs. (26, 27) can be considered as linear homo- 
geneous algebraic equations for s\n.{v + (3) and cos{v + (3) . 
They have a solution if 



d'9 d'& 
dui du2 



+ 



dip dip 



cos^ ip dui du2 



0. 



(33) 



But if we make an assumption given by Eq. (31) then 
Eq. (33) results in 



dp dp 
dui du2 



COS^ ip 



dp 
dp 



= 0. 



(34) 



If Eq. (32) holds, it follows from Eq. (34) that p = 
const and § = const, i.e., there is no variation of the 
electromagnetic fields in space and time. 

The fact that wc have not found other pulse group 
velocities than c and (c^^ -|- G/fi-^) is in full agreement 
with the results of perturbative approach [16,17]. 



IV. NUMERICAL SOLUTIONS 

It is natural to expect that any pulse of finite dura- 
tion evolves in the medium into pair of fast and slow 
pulses, which become more and more separated in space 
due to the difference of their group velocities. Indeed, 
our numerical simulations confirm such an expectation. 
An example is shown in Fig. 2. The quantity wq used for 
normalization of the horizontal axes of the plots in Fig. 2 
and subsequent determines the order of magnitude of ip 
and d, which are ~ /{Gwq). The adiabaticity condi- 
tion (8) results in the following rcstric;tion: wq ^ Q./G. 
One can see that the incident pulse evolves into a well 
separated pair of fast (F) and slow (S) pulses, and the 
mixing angle v describing transitions between the two 
dark states emerges (the incident pulse is chosen in such 
a form that i/ = at the medium entrance) . The F and S 
pulses at large propagation distances can be excellently 
fitted with formulae (30) and (29), respectively. 

We also investigated numerically collisions between 
fast and slow pulses. The results are presented in Fig. 3. 
The pulse sequence is organized in such a way that the 
pulse of a shape satisfying Eq. (29) enters the medium 
first. After some time delay the next pulse obeying 
Eq. (30) enters the medium. The first pulses propagates 
at the slow group velocity whereas the second one prop- 
agates at the speed of light. The distance between them 
decreases, and at certain Q the two pulses overlap (this 
is marked by O in Fig. 3b). Their nonlinear interaction 
leads to strong distortion of their shapes, which becomes 
apparent at larger propagation distances. Thus adiabatic 
pulses in a tripod medium cannot be called solitons in the 
exact sense of soliton definition by Zabusky and Kruskal 
[22]. Note that it is impossible to arrange a collision of 
two adiabatons in a A-medium. 

The parameter /3 is equal to 1.12 for Fig. 2 and 1.87 
for Fig. 3. 
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FIG. 2. Splitting of the incident pulse into the fast (F) and 
slow (S) pulses, (a) Sine functions of the angles a = i? (thick 
line) and ip (thin line) versus scaled nonlinear time, w/wo, at 
the medium entrance, C = (the boundary conditions); the 
same line styles are reserved for i? and (p in all the subsequent 
figures, (b) The same as in (a) for — 20wo (numerical 
simulation results), (c) Sine of mixing angle v versus w/wo 
for C = (dashed line) and 20wo (dotted line) . (d) Fitting of 
the plot (b) with Eq. (30) (solid lines; Ci = 0.88, C2 = 1.75) 
and Eq. (29) (dashed lines; C3 = 0.43, d = -0.90. Axes are 
in dimensionless units in all plots. 
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FIG. 3. Collision of pulses, (a) Sine functions of the angles 
a — & and ifi versus w/wo at = (the boundary conditions), 
(b) The same as in (a) for ^ = lOwo (numerical simulation 
results), (c) The same as in (b) for ^ = 30wo. See the text 
for more detailed description. 



V. PARTICULAR REGIMES OF PROPAGATION 

There are a few particular regimes of adiabatic pulse 
propagation admitting analytic treatment. The first 
one occurs if atoms are prepared initially in a sta- 
tistical mixture described by the density matrix g = 
5 + rather than in a pure 

state. Such a mixed state remains invariant under the 
action of action of the slowly varying electromagnetic 
fields: BqW = g, where B is given by Eq. (18). Sta- 
tistical averaging over g is equivalent to averaging over 
the parameter (3 uniformly distributed between and 27r, 
without any correlation with the instantaneous values of 

and 1?. The result of statistical averaging of Eqs. (22) 

is 



dip 



2'duj' 



1 d%') 

2 dw 



(35) 



Equations for ip and ^ become decoupled. Their solu- 
tion = ipo{'w — I (J), 1? = ■!?o(w — i() describes inde- 
pendent propagation of perturbations of and ip at the 
same group velocity Vg = \c~^ + G/(2ri^)] 
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FIG. 4. Adiabatic pulse propagation for the particular ini- 
tial condition (3 = 0. (a) Sine functions of the angles a = •& 
and ip versus w/wo at C = (the boundary conditions), (b) 
The same as in (a) for ^ = 20wo (numerical simulation re- 
sults) . 

We may hazard a conjecture what occurs if atoms are 
prepared in a mixed state with the density matrix g' = 
g|$(i))($(i)| +(i_g)|$(2)>^$(2)|^ < g < 1. It is 

Ukely that there are always two classes of pulses with 
well defined group velocities. If q grows from to 0.5, one 
of these velocities decreases whereas the other increases. 
At q = 0.5 they achieve the same value mentioned in the 
previous paragraph, and then again restore their values 
c and (c~^ + G/n'^)~^ , as q approaches 1. At least, it 
can be proven easily in the pcrturbative regime, when 
the changes of both ip and 'd arc small. 

Another interesting regime is related to particular ini- 
tial conditions /3 = or /3 = 7r/2. Let all the atoms 
be pumped initially into the state |3). The fields are 
switched on in the following order, which is a generaliza- 
tion of the counterintuitive pulse order for a A-medium 
[3]: Initially, at w = only the field driving the empty 
transition |2) |0) is present, i.e., = = 0. Obviously, 
/? = 0. Then the field driving the transition |1) •(-> |0) is 
switched on adiabatically, so that ?? grows and then is 
kept constant at a certain level. Finally, the field driving 
the transition |3) |0) is switched on. 

When 1? changes, simp = 0. Then, according to 
Eq. (17), 1/ remains zero, and Eqs. (25) are reduced to 
d(p/ {dui) = 0, di)/ {du2) = 0. Such a propagation regime 
occurs unless the front of the i?-pulse, propagating at the 
slow group velocity, approaches the front of the (/3-pulse, 
propagating at c. 

Thus one has a possibility of preparation of a tripod 
medium in any desired coherent superposition of low- 
energy states. Numerical results presented in Fig. 3 il- 
lustrate this conclusion: Finally, atoms in the region 
< ^; < 20wo are prepared in the state —0.29 |1) — 



0.53 |2) — 0.80 |3), as can be derived from the values 
= —0.65, ■& = 0.50 'At w = AOwq. Then one can sud- 
denly change the laser radiation parameters in such a 
manner that this state will correspond to a coherent su- 
perposition of the two dark states defined with respect 
to the new values of the Rabi frequencies, thus obtaining 
a new value for the parameter (i. 

The case of /? = 7r/2 is physically equivalent to the 
previous one, differing only in notation of the states and 
electromagnetic fields. 



VI. CONCLUSION 

Requirements for experimental implementation of adi- 
abatic pulse propagation in a tripod medium should 
not differ from that for slow light propagation in A- 
media [2,4,8]. A method for initial preparation of a 
tripod medium in any desired superposition state was 
outlined in the previous section. For example, con- 
sider a tripod medium with the following parameters: 
fioj ~ 10~^* esu • cm, kj k. 10^ cm~^, n « 10^^ cm~^. Let 
the total laser intensity be of about 3 mW/cm^ (slightly 
below the atomic transition saturation limit). Hence, 

« 2.4 • 10^ s-i and Vfi/G » 10* cm/s (i.e., the group 
velocity of the slow pulses is « 0.3 • 10~^c). The value 
of the scaling parameter of the horizontal axes of the 
Figs. 2 - 4 uiQ « 0.1 cm is thus large enough to pro- 
vide adiabaticity. Therefore the processes illustrated in 
Fig. 2-4 can be observed in a few centimeter long gas 
cell. The time delay between the fast and slow pulses is 
of about 0.3 ms, therefore the lifetime of coherence be- 
tween the states | j), j = 1,2,3, should be 1 ms or longer. 
It is achievable in coated cells or cells with a buffer gas. 

To conclude, we have investigated electromagnetic 
pulse propagation in a coherent atomic medium with the 
tripod configuration of levels in the adiabatic regime. 
The propagation equations (26 - 28) are derived and 
their solutions in the form of slow [Eq. (29)] and fast 
[Eq. (30)] pulses are obtained analytically. Our numer- 
ical simulations confirm that these solutions are general 
asymptotic solutions for any incident pulse of a finite du- 
ration. We have suggested a method of preparation of a 
tripod medium in an arbitrary superposition of the low- 
energy states based on switching on the laser fields in a 
counterintuitive order. The tripod scheme provides two 
novel features in comparison to the A-scheme. The first 
one is adiabatic pulse propagation in a medium prepared 
in a statistical mixture of the two dark states. The sec- 
ond one is the possibility of collisions between the slow 
and fast pulses revealing that they change their shapes 
after nonlinear interaction and thus do not satisfy the 
classical definition of a soliton [22] . 
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